Observing the temperature of the Big Bang through large scale structure 
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It is widely accepted that the Universe underwent a period of thermal equilibrium at very early 
times. One expects a residue of this primordial state to be imprinted on the large scale structure of 
space time. In this paper we study the morphology of this thermal residue in a universe whose early 
dynamics is governed by a scalar field. We calculate the amplitude of fluctuations on large scales 
and compare it to the imprint of vacuum fluctuations. We then use the observed power spectrum of 
fluctuations on the cosmic microwave background to place a constraint on the temperature of the 
Universe before and during inflation. We also present an alternative scenario where the fluctuations 
are predominantly thermal and near scale-invariant. 

PACS numbers: 0000000 



A cornerstone of modern cosmology is that the universe 
underwent a sustained period of thermal equilibrium at 
early times. Two of the key predictions of the big bang 
cosmology, the spectrum of the cosmic microwave back- 
ground and the abundance of light elements, hinge on the 
existence of this primordial hot phase [l[ . A key charac- 
teristic of systems in thermal equilibrium is the presence 
of fluctuations. These are, to some extent, uniquely de- 
fined and can be derived from the microphysical prop- 
erties of the system @- Hence we expect to be able to 
characterize the fluctuations of the energy density of the 
early universe which in turn lead to irregularities in the 
fabric in space time. These should be reflected in the 
distribution of large scale structure, the propagation of 
light rays and other such cosmological observables. 

There have been a number of attempts at pinning down 
the fine details of the thermal fluctuations in the early 
universe. Under standard assumptions it can be shown 
that thermal models are observationally unsound. A 
generic n s = 4 prediction for the spectral index follows, 
unless there is a phase transition, in which case n s = 0. 
This can be bypassed, and a more congenial n s w 1 be 
predicted, by considering non-standard assumptions: e.g. 
by considering a gas of strings at the Hagedorn phase [3| , 
or by invoking an early holographic phase in loop quan- 
tum cosmology Q, followed by a phase transition. One 
can also appeal to the technicalities of loop quantum cos- 
mology 5] or postulate a mildly sub-extensive contribu- 
tion to the energy density Q • All these scenarios require 
speculative new physics. 

In this paper, we revisit this issue by focusing on what 
has become a standard and fruitful model of the universe: 
a perturbed homogeneous and isotropic spacetime whose 
dynamics is driven by a scalar field. Without loss of gen- 
erality, we will restrict ourselves to a scalar field with 
an exponential potential but will allow both positive and 
negative kinetic energies . If the field rolls sufficiently 
slowly away from the origin, we have power law, acceler- 
ated expansion. If the field rolls sufficiently quickly, the 



energy density in the scalar field will mimic the behaviour 
of an assortment of cosmological fluids (such as radiation 
or dust) . If the kinetic energy of the scalar field is neg- 
ative, we obtain "phantom" -like behaviour: the effective 
equation of state w = P j p (where p and P are the en- 
ergy density and pressure in the scalar field) is such that 
w < — 1. Such a setup allows us to analytically calculate 
the amplitude and spectrum of thermal fluctuations in- 
cluding gravitational backreaction. In this paper we will 
focus on universes that underwent superluminal expan- 
sion. 

Let us briefly revisit the model. We will consider 
a potential for the scalar field of the form: V(4>) — 

exp(— y^j^j) where Mpj is the reduced Planck 

mass. The evolution of the scal ar field is given by 
4> = \f2pMpi Ya{Mpxt/ \/p(3p — 1)) and the Friedman 
equations lead to a simple solution of the form a oc t p 
and H = da — p/t where ■ = d/dt. Note that, if p > 1, 
the expansion is superluminal. It is convenient to rewrite 
some of these results in terms of conformal time, r. If 
p > 1 we have that the past is at r = — oo and blows up at 
t = 0. We then have that the scale factor and the confor- 
mal Hubble parameter is given by a oc (— t)"- 1 and the 
conformal Hubble parameter is given by 7i = - 
where ' = djdr. 

Let us now focus on how perturbations on these back- 
ground cosmologies are seeded and evolve [Io|. Recall 
that we can expand a scalar field and space-time met- 
ric around a homogeneous background, <j> = cf>o + tp and 
ds 2 = a 2 [(l + 2$)dr 2 - (1 - 2*)dr 2 ]. The quantity of 
choice is the gauge invariant variable, 

v = a(5<p + ^) 
ti 

which can be related to the curvature perturbation, TZ = 

— v / z where z = ^W- The gauge invariant Newtonian (or 
"Bardeen") potential, $, can be found from 1Z through 
fc 2 $ = iwGcbozlZ' For our choice of background cosmolo- 
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gies, the gauge invariant perturbation variable obeys a 
Besscl equation with a general solution given by: 

Mr) = M\r\) 1/2 Ju{k\r\) + B^r^Y^k^) 

where J v (x) and Y v (x) are Bessel functions with v = 
— -i- — — — 

2 ' p-1 • 

In a universe undergoing superluminal expansion there 
is a natural mechanism by which fluctuations can be 
seeded. We assume that v is promoted to a quantum 
operator: 

f <i 3 k 

u(t,x)= J j2^[v k (T)a k e*-* + vt(T)ale- ik -*} 

where Sk and its conjugate are the annihilation and cre- 
ation operators for the n particle state. We are interested 
in the two point correlation function (£(x+r)?)(x)), where 
(• • •) is a quantum expectation value depends on the state 
one is considering, i.e. (A) = A natural choice is 

the ground-state or vacuum state of each mode, \Q = |0). 
In the past, where (—kr) — ► — oo, a given mode was well 
within the horizon. This allows us to uniquely define the 
solution (i.e. the coefficients and £>k) to the mode 
equation to be v k (r) = ^ e ^+ 1 / 2 ^l 2 {-T) 1 ' 2 H^ ] \-kr) 

(where H„(x) is a Hankel function). This solution has 
a unique behaviour at late times (i.e. when (— kr) — ► 0): 

1H M -> P *("-l/2)ir/2 2 i'-3/2 F M 1 f i \-i/+l/2 

Wk(T) e 2 r[3/2]V2fc ( fcT) 

1 4(1 -u) ] 

1Z can be trivially obtained from the above solution and 
we find that 

2vr2 ^r[3/2]J Z 2 fc I 1 ) 

which goes to a constant as p — * oo. We find the well 
known result that the scalar spectral index is given by 
jiq — 1 = -r 2 — = ^fT7> ■ We can see that, in the limit of 
to -1 we have pure scale invariance. 

Throughout the above calculation we have discarded 
any reference to the hot origins of the universe. Yet we 
are starting off at high energies, when the Universe would 
have been strongly interacting. It would be natural to ex- 
pect the imprint of these thermal initial conditions on the 
scalar field in some way. Indeed one would expect fluc- 
tuations in the scalar field to be thermalized through a 
variety of different mechanisms. The universe may have 
entered a scalar field dominated regime from a preceding 
radiation dominated regime; interactions with the hot ra- 
diation would have led the fluctuations in the scalar field 
to be thermal. Furthermore, the scalar field model we 
are considering has non-linear self-interactions through 
the exponential potential. Very short wave modes would 
play the role of a heat bath even through the period of 



superluminal expansion and scalar field domination. The 
details of how primordial fields undergo evolution in a hot 
phase have been studied in great detail in [13, EH where 
a number of effects where identified emerging from the 
non-equilibrium nature of the problem. 

In what follows, we will disregard non-equilibrium 
effects: these will introduce small corrections and 
can be included in a more detail calculation. Our 
calculation is therefore undertaken in the setting of 
equilibrium statistical mechanics: the appropriate ex- 
pectation value to consider is given by (A) = 
Y, n Pnn{n\A\n)/(Y Jn Pnn{n\n)), where \n) is the n- 
particle state (referring to a given momentum k). The 
simplest approach is to simply posit that each mode is 
Boltzman weighted. Recall that this involves setting 
the density matrix above to p nn — e~@ En where E n is 
the energy of a given mode with occupation number n, 
j3 = 1/KbT, Kb is the Boltzman constant and T is the 
temperature. Hence we find that 

^-^\ Vk (r)\ 2 [2n(k) + l]e^ r 

where the resulting number density (subtracting out the 
vacuum state) for is given by n(k) — e)3B(fc 1 T) _ 1 . 

The energy of the perturbation can be found from 
the Hamiltonian density of v. With the above solu- 
tions we have that E(k,r) = h7T ^ T \ B(k, r)\H v {k\T\) with 
B(k,r) = (3p 2 -p)/[t 2 (p~ l) 2 ] + k 2 . Note that in the 
short wave length limit we recover the standard for plane 
waves: E(k,r) — hk. Also note that the energy is de- 
fined in terms of conformal quantities (derivatives are 
taken with regards to r and conformal x). This means 
that the temperature we use must also be in the same 
conformal frame: T — T p h ys a where T p h ys is the physical 
temperature. If we assume that thermalization is main- 
tained through a heat bath which evolves as radiation, 
and there is no generation of entropy during inflation, we 
have that T = constant. This value of T will be crucial 
in what follows. 

We have not yet arrived at the final result. Any 
mechanism that keeps the scalar field in thermal equilib- 
rium must break down as a given mode becomes larger 
than any causal scale, i.e. around the horizon scale. 
On superhorizon scales, we expect the spectrum to be 
frozen in — the heat bath or interactions are irrelevant. 
In other words n(k) will be frozen at the value it has 
when k\r\ ~ 1. This means that the above expression is 
not entirely accurate and we must replace E{k, r) with 
T")|feT|~i- We can now express our main result: the 
full spectrum of fluctuations, including the thermal con- 
tribution is 

protal ik) =v Th {k)+v L{k) = pO (fc)[2(n(fc>T) +1] 

_ (2) 

This result is slightly different from that in |6[ (the model 
is not the same). It also cannot be directly compared 
with the results of warm inflation [l3| . 
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Let us now explore the consequences of Equation [2] 
As a first guess, one would expect to be in the Rayleigh- 
Jeans regime when the modes exit the horizon. We then 
have, on superhorizon scales, 



(p-l) 2 (-T)-^ +1 k^ 



Ap 2 - 3p + 1 
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If we reexpress Eq. [3] in terms of the equation of state, 
we have that 



n. , - 1 = 



5 + 3w 
1 + 3w 



Close to de-Sitter we find that n s ~ 0, that is white noise: 
we do not get scale-invariance because the temperature 
is decreasing like 1/a, breaking the deSitter invariance 
(this is to be contrasted with the work of [ll|). Instead 
we find that a scale invariant spectrum arises if we assume 
a "phantom" regime with w — —5/3. 

Our expression is insensitive to the details of thermal- 
ization and horizon crossing and it gives us a reasonable 
idea of what to expect. A useful exercise is to compare 
the contribution of thermal fluctuations relative to vac- 
uum fluctuations during an inflationary period. We have 
that 



1 n 

' K 



(p - l) 2 



K B T 



4p 2 

For p> 1 we find 
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In general the prediction of this model is a break in the 
power spectrum at pivot scale k p ~ 0.1T. For k < k p the 
fluctuations are predominantly thermal with spectral in- 
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- 1, with the quantum fluctuations spectral 



index n° s given by the usual formula. In this regime we 
are invariably in the Rayleigh Jeans limit. For k > k p 
the fluctuations are predominantly quantum, with ther- 
mal fluctuations suppressed by a factor of e~ k / kp , given 
that we are in the Wien regime. At horizon crossing we 
always have E ~ k, so this can be replaced in the formula 
for n{k) in either regime. 

We now examine the implications of this result for 
two viable scenarios, where the fluctuations are predomi- 
nantly quantum and thermal, respectively. If we have an 
inflationary scenario (w « — 1) then the dominant fluctu- 
ations on observable scales should be quantum, for these 
are near-scale-invariant. A priori the prediction of this 
model is a turn over in the spectral index from n s 1 to 
n s = on large scales (for k < k p ). There is clearly no 
evidence for higher power in the lowest multipoles of the 
CMB so, at best k p could be the current horizon scale 
kno- This is reflected on an upper bound on the tem- 
perature during and before inflation or alternatively on 
a constraint on the ratio of the temperature before and 
after reheating (T& and T a ). Recall that the conformal 



temperature T = T p h ys o. is only a constant if there is no 
entropy production, so that it does suffer a jump, from 
T to T a at reheating. Bearing this in mind, k p ~ 0.1X& 
but kn ~ aH(). Therefore k p < kno translates into 
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(4) 



A marginally tighter bound can probably be obtained 
through the Grischuk-Zeldovich effect: superhorizon fluc- 
tuations with such a red spectrum will further boost the 
quadropole (TH- We can convert our constraint into a 
physical temperature during inflation if we assume a spe- 
cific model. For example, if the inflation ended at the 
GUT scale, when the energy scale is of order 10 17 GeV 
and at a redshift of z ~ 10 28 , the temperature of the 
Universe just before reheating would have been, at most, 
10 -2 eV. This means that the Universe hits the Planck 
temperature more than 68 efoldings before reheating, so 
that there is scope for producing the observed structure 
of the Universe (for which 50 to 60 efoldings before re- 
heating is enough), but, if the bound is saturated, not 
much more. In general the bound Q forces the maxi- 
mum number of efoldings to be 
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Epi 

Elnf 



2.3 



(5) 



If we can assume that H doesn't vary by much during 
inflation, and if all the energy in the inflaton field is con- 
verted into radiation during reheating we can translate 
the bound ((4]) into M > 64. Relaxing these assumptions 
produces a tighter bound. This seems to rule out open 
inflationary models. 

If we have a phantom scenario with w ~ —5/3, the ob- 
served structure of the Universe should be thermal. The 
prediction is a near scale invariant spectrum breaking 
into n s = 2 for k > k p . Thus we should have k p > kso, 
where kso is the smallest scale for which the primordial 
power spectrum is observable. The constraint is now an 
upper bound on how much entropy has been produced 
since the observed structure left the horizon; specifically: 



T b _ hk S o 
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(6) 



where we have assumed that the smallest scales that can 
be probed are of the order of a Kpc. In this scenario 
we have roughly that a oc l/(— t), with t < (i.e. p — 
— 1), p oc a 2 , and H^ 1 = —t (the horizon's physical size 
decreases) . Also a = — 1/ (2Ht) . 

The normalization in this model is obtained from a 
constant of motion combining the energy in the ther- 
mal bath and that in the background field. The rele- 
vant factor in ([3]) is 1/(z 2 t) which can be rewritten into 



T phys H/M%i, i.e. 



1/4 1/2 



or - T phys /\t\. It's this im- 



portant constant that must be ~ 10~ 10 to match obser- 
vations. Should all the energy in the "phantom" field be 
converted into radiation at the end of this phase we there- 
fore get the rather undemanding bound T a < \Q A Mpi (in 
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combination with ([6])). But by requiring that the cur- 
rent Hubble volume was once inside the phantom Hub- 
ble volume (in a calculation mimicking the inflationary 
counterpart) we find that T phys /T Pl ~ 0.1 at the start of 
the phantom phase (and that requires saturating bound 
©). If the thermal bath is set up at TP^ S - T Pi the 
break into n s = 2 should happen only an order of mag- 
nitude or so above kso- Whether this could be observed 
is debatable. 

Note that for simplicity we have considered w — —5/3, 
but strict scale-invariance in this scenario is actually 
pathological as it requires it is only for— 5/3 < w < —1/3 
that the Newtonian potential <& stays constant and has 
the same spectrum as 1Z on large scales. For w < —5/3 
the potential diverges. However as long as the spectrum 
is slightly red this is not a problem and we have for the 
growing mode: 



We conclude with a few comments on aspects of this 
model, and how they relate to other work. We stress that 
our system is very different from a single thermal fluid, 
as previously studied 0, H, ■ Here the unperturbed field 
<fro is not thermalized; only its fluctuations 5<f> are ther- 
malized. The fluid </>o drives the expansion and provides 
the leading order energy, but no entropy. Whatever the 
equation of state w for <fio, the 5<f) behave like standard 
thermal radiation, with w± = 1/3 and supply the en- 
tirety of the entropy of the system. This feature allows 
us to bypass a number of thermodynamical constraints 
pertaining to single thermal fluids, namely the relation 
( — 1 + 1/w between the ( exponent appearing in p cx 
and w. If we insist on a Stephan-Boltzman law of the 
form pq cx (where po is the energy in 0o, and T is the 
temperature of 8<p) we find instead that ( = 3(1 + w). 
This doesn't contradict any fundamental thermodynam- 
ical constraint: the usual result merely indicates that the 
second order energy, contained in 8(f), should go like T 4 . 

But even a two-fluid model breaks down when dis- 
cussing thermal fluctuations. Indeed Maxwell's formula, 
c%(R) — T 2 dUJdT, which is the workhorse of much pre- 
vious work [3, 4, 5], is not applicable here. The energy 
fluctuation is of the form Sp ~ (f>o8<p, i.e. a cross term be- 
tween the unthermalized 0o and the thermalized 5cj). So 
the energy fluctuation of the system is, to leading order, 



<y%{R) cx UqUi (-R), where Uq = poV is the average energy 
in (fig, and U\(R) is the average energy in 5<fi smoothed on 
scale R (which is ~ T). Unusually, we only need to know 
the average energy of the thermalized system to work out 
the leading order energy fluctuation in the overall system. 
These novelties conjure to bypass the general prediction 
n s = 4, allowing for scale-invariant thermal fluctuations 
without appealing to any new physics. 

Regarding the Gaussianity of these fluctuations it has 
been shown Q that for a single thermal fluid ther- 
mal fluctuations are very approximately Gaussian in the 
Ray leigh- Jeans limit (but not in the Wien limit). How- 
ever, just as it happens with the equivalent calculation of 
the variance, the calculation of the cumulants in a single 
thermal fluid is not applicable to our system. Instead 
we note that the derivation of Gaussianity usually used 
for linear inflation applies to any density matrix that is 
diagonal in the number operator, including a thermal 
state. We therefore expect the thermal component to be 
Gaussian, too, rendering the thermal scenario presented 
above viable. This is in contrast with non-linear infla- 
tionary couplings, that may produce a certain degree of 
non-Gaussianity [151 ]. 

Finally, we remind the reader that we are considering 
a universe that starts off in thermal equilibrium. The 
hallowed example is that of what has become known as 
new Inflation: as the Universe cools down, the scalar field 
settles down into a slow roll regime and it is potential 
energy dominated. This is not, however, a generic feature 
of the inflationary cosmology. One appealing alternative 
is a Universe that emer ges through quantum tunnelling 
into an inflationary era [16]. Another possibility is that 
our local patch has entered into an inflationary regime 
as a result of a Planck scale fluctuation of the Inflaton 
[l7j . The initial state for the onset inflation would not 
necessarily be thermal. In both of these scenarios we 
don't expect a thermal imprint on space time on large 
scales. 
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